Nernst effect in the two-dimensional XY model 
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We calculate the Nernst signal directly in the phenomenological two-dimensional XY model. 
The obtained numerical results are consistent with the experimental observations in some high-Tc 
cuprate superconductors qualitatively, where the vortex Nernst signal has a characteristic tilt-hill 
profile. It is suggested that the excitations of vortex and anti- vortex in 2D is the possible origin of 
the anomalous Nernst effect. 

PACS numbers: 74.25.Fy, 74.40.+k,74.72.-h 



By a phase only model, Podolsky et al[lj have studied 
the Nernst effect [21 and diamagnetic response due to the 
thermal superconducting fluctuation through the trans- 
verse thermoelectric conductivity (axy)- They found that 
the Nernst effect has a much sharper temperature de- 
cay than predicted by Gaussian fluctuations. We think 
that for this phase model the Nernst signal, cn = ctxyjo 
{a is the electrical conductivity), can be calculated di- 
rectly and can be compared directly with the experi- 
ments. For some cuprate superconductors, such as single- 
layer Bi2Sr2-yLayCuOe, the quasi-particle contribution 
to the Nernst signal is virtually negligible for y in the 
range 0.4 — 0.6. The phase model can only be ap- 
plied to such kind of cuprate superconductors. Note that 
c^xy = cXxy + '^xy^ which includes the contributions from 
both vortices and quasi-particle, and can not be mea- 
sured directly in experiments. It is also difficult to sepa- 
rate the contributions from vortices and quasi-particles. 
The contribution from quasi-particles a^y can not be ob- 
tained in this model. Therefore the axy in Ref. [l| in 
principle can only describe the vortex dominated Nernst 
effect, and can not describe the Nernst effect in those 
cuprate superconductors where the quasi-particle signal 
is comparable with the vortex one. 

The same Hamiltonian of the two dimensional (2D) XY 
model is given bytl] H = ~JY,{ij) cos(0i - (jjj - A^j), 
where Aij represents the magnetic vector potential of a 
field B = V X A perpendicular to xy plane. The average 
number of vortex lines per plaquette / = a^B/^Q, with a 
grid spacing. With the model-A dynamics, the equations 
of motion are given by — ~^§j"^^ii't) where T is 

a dimensionless constant which determines the time scale 
of relaxation. The thermal noise term Ti (t) is assumed 
to satisfy {r^{t)) = and {r^{t)Tj{0)) = 2hTkBT5^■j5{t). 
The time and the temperature are in units of h/TJ and 
J/ks- The present simulations are performed with the 
system size L = 64 for two directions. 

The temperature gradient is set up along x direction, 
where a open boundary condition should be applied. 
Specifically, we introduces a global "twist" variable Aj,(t) 
to track the average phase drop per link in the y-direction 
[3, Q • The voltage by the diffusion of free vortices ther- 
mally driven is given by Vy = Aj^. Since the periodic 



boundary condition for the phase variable is satisfied in 
the y -direction, previous pseudo-spectral algorithm Q 
can also be employed to facilitate the simulations. The 
time stepping is done using a second-order Runge-Kutta 
scheme with At = 0.05. The measurement is performed 
after the steady-state is reached. It is found that the 
present thermally driven system is time-consuming con- 
siderably. Our runs are typically (4 — 8) x 10^ time steps. 
The detailed procedure in the simulations will be de- 
scribed in a full paper @ . 

We calculate the Nernst signal directly through ejv = 
1^ "^1 . As shown in Fig. 1(a), our numerical results are 
consistent with those in Figs. 11 and 12 in Ref. Q 
qualitatively, where the vortex Nernst signal has a char- 
acteristic tilt-hill profile. At each temperature, in the 
low magnetic field regime, the interaction between vor- 
tices is small and can be neglected, so the Nernst sig- 
nal is added linearly with the addition of free vortices 
as the fields increase. In the high magnetic field regime, 
the strong interaction among the dense vortices prohibits 
the vortex motion driven thermally, and cause very small 
Nernst signal. The competition of these two mechanisms 
results in a maximum Nernst signal in the intermediate 
fields. This nonmonotonous behavior is clearly observed 
in the curves for the Nernst signal at various tempera- 
tures above and below Tkt in the present simulations. 
For fixed fields, as the temperature increases, the Nernst 
signal also displays nonmonotonous behavior, also sim- 
ilar to the experimental observations [1| . It is suggested 
that the excitations of vortex and anti- vortex in 2D is the 
possible origin of the anomalous Nernst effect observed in 
strongly layered high-Tc cuprate superconductors. The 
inter-layer Josephson coupling only renormalizes the in- 
plane coupling strength and may play no essential role in 
this phase-disordering scenario. 

More interestingly, our numerical results for the Nernst 
signal can explain the ridge flelds in the Nernst experi- 
ments well [3] . In this model, the critical field Hc2 is cor- 
responding to / = 1, i.e. the vortices occupy the system 
completely. The peak position around / = 0.1 shown in 
Fig. 1(a) is consistent with the typical ridge fields Q.\Hc2 
observed in the Nernst experiments 3 • 

In Fig. 1 of Ref. [l|, one can find that the transverse 
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FIG. 1: (a) CAT and (b) function of / at various 

temperatures above and below Tkt = 0.89. The system size 
L = 64. 



thermoelectric conductivity a^y decreases monotonically 
with temperatures. In our simulations, we find a non- 
monotonous behavior of a^y below Tkt, as shown in 
Fig. 1(b). This numerical prediction for the only con- 
tribution from vortices in axy is clearly called for further 
confirmation in further advanced experiments where the 
vortex contribution is indeed separated from the total 
Uxy To the best of our knowledge, at the moment no 
one can separate the vortex contribution easily. 



In short, the direct simulation for the Nernst signal in 
the phenomenological 2D XY model can provide a quali- 
tatively explanation to the experimental findings in some 
high-Tc cuprate superconductors. 
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